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We investigate the electrostatics of a partially formed, charged spherical shell in a salt solution. 
We solve the problem numerically at the Poisson-Boltzmann level and analytically in the Debye- 
Hiickel regime. From the results on energetics of partially formed shells we examine the stability of 
tethered (crystalline) and fluid shells towards rupture. We delineate different regimes of stability, 
where, for fluid shells, we also include the effects of bending elasticity of the shells. Our analysis 
shows how charging of the shell induces its instability towards rupture but also provides insight 
regarding growth of charged shells. 
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I. INTRODUCTION 

Self-assembly as well as disassembly of biological 
macromolecules and their aggregates is governed by the 
standard principles of colloid and nanoscale stability the- 
ory [l| : repulsive interactions of in general disparate ori- 
gins are usually counteracted by the ubiquitous van der 
Waals (vdW) force 0. For macromolecular aggregates 
bearing net charges, the Coulomb interaction in one of 
its many guises [3|, [j] is the leading component of repul- 
sive interactions, while for those decorated with exposed 
hydrophobic moieties the hydrophobic interaction [5(, at 
least in part related to the vdW interactions, joins in the 
overall attractive part of the stability condition. A stable 
state of (classical) matter is possible only if other forces 
apart from electrostatic ones are present Q . In addition 
it was shown recently [7[ that similarly the vdW interac- 
tions alone, at least for bodies with permittivities that are 
all higher or lower than that of the intervening medium, 
are also not able to generate stable states of matter. Only 
by combining the two together, or more unusually having 
bodies with alternating higher/lower permeabilities than 
the intervening medium, can full stability of matter be 
achieved. 

In biological context it is mostly the screened electro- 
static interactions and hydrophobic/van der Waals in- 
teractions that set the structural stability conditions of 
macromolecular aggregates Of the two mentioned 
above, the hydrophobic force corresponds at least in part 
to an effective interaction arising from the statistical 
properties of water molecules around the hydrophobic 
moieties @] , while the van der Waals force is a true force 
which can be modified by the presence of the aqueous sol- 
vent but acts also in its absence. If the attractive forces 
are strong enough to overcome the electrostatic repul- 
sion, the macromolecules that in general contain charged 
as well as hydrophobic patches may cluster into different 
stable aggregate structures, depending on their geome- 
try as well as on the spatial distribution of hydrophobic 



and hydrophilic patches along their surface. Since the 
presence of an aqueous solution, which in the biologi- 
cal milieu contains various dissolved ionic species, affects 
mostly the properties of the electrostatic interaction via 
a pronounced Debye-Hiickel screening, the process of ag- 
gregation itself is also influenced by the screening prop- 
erties of the solvent, i.e. the concentration and valency 
of salt ions. 

An interesting example of a soft matter system where 
the electrostatic interactions determine its stability is the 
problem of the shape of a charged fluid droplet studied by 
Lord Rayleigh more than a century ago Q. He showed 
that a charged droplet has a limited range of stability 
and breaks into smaller droplets if its charge exceeds a 
certain critical value. This critical value of the charge is 
set by the balance between its electrostatic energy, pro- 
portional to the square of the charge, and its surface en- 
ergy, proportional to the surface tension of the air-water 
interface. The latter can be associated with hydropho- 
bic interactions at the fluid-gas interface. Beyond the 
limit of stability set by the critical value of the charge, 
the equilibrium configuration is a collection of smaller 
droplets with charges below the critical value, placed at 
an infinite distance from each other. 

The concept of Rayleigh's instability in charged 
droplets described above can also be applied to polyelec- 
trolyte globules in poor solvents, which again implies the 
existence of hydrophobic attractions between monomers. 
The stability limit in this case depends on the fraction of 
charged monomers in the polyelectrolyte chain. Just as 
in the case of a charged droplet, a polyelectrolyte chain 
in a poor solvent can reduce its energy by splitting into 
a set of smaller charged globules. Because of the connec- 
tivity of the polymer chain, the smaller globules have to 
remain connected by an intervening finite polymer seg- 
ments (To| . In such a way the equilibrium structure after 
the limit of stability of the original polymer globule is 
similar to a pearl necklace. Beyond the limit of stability 
the polyelectrolyte chain first takes a dumbbell configura- 
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tion, turning into a pearl necklace with three beads joined 
by two segments, leading eventually to a whole cascade 
of transitions between necklaces with different number of 
beads as the charge on the chain increases fill ]. 



The system that we study in this paper belongs to 
the class of systems discussed above. In what follows, 
we will be interested in stability of spherical, shell- 
like macromolecular aggregates composed of charged 
(macro)molecules. What we have specifically in mind 
are two rather dissimilar systems: on the one hand soft 
spheroidal charged lipid vesicles and on the other rela- 
tively hard icosahedral empty virus shells (capsids). In 
the former case the shell is composed of a quasi 2D fluid 
layer of lipid molecules [12j . whereas in the latter the 
shell can be envisioned as a crystal-like assembly of teth- 
ered capsomeres [l]|. Once assembled, such structures 
may become unstable if the amount of salt dissolved in 
the bathing solution, quantified by its ionic strength, is 
decreased, so that the screening of electrostatic repulsion 
becomes less efficient flij . For the macromolecular ag- 
gregates composed of complex molecules such as proteins, 
as in the case of virus shells, the amount of charge they 
carry may be modified by shifting the dissociation equi- 
librium via the solution pH [l5[ . Apart from the changes 
in the ionic strength this presents yet another effective 
way to influence the stability of shell-like aggregates. In 
general, once the shell becomes unstable due to changes 
in the electrostatic interactions brought on by either the 
changes in ionic strength or pH, it may break up accord- 
ing to different scenarios. We shall consider here only the 
conceptually simplest cases, and we shall show how the 
consideration of these simple cases leads to more general 
conclusions regarding shell stability. We shall apply our 
results both to fluid and tethered charged shells, so our 
findings should be of importance in the study of stability 
of charged fluid vesicles as well as virus capsids. 



The structure of the paper is as follows: we shall first 
introduce a simplified model of partially formed fluid and 
tethered spherical shells (Sec. Ill A|) . We shall then for- 
mulate the theory of electrostatic interactions for such 
simplified geometries based on the solution of either the 
complete Poisson-Boltzmann (PB) theory (Sec. Ill Bp or 
its simplified derivative in the Debye-Hiickel (DH) form 
(Sec. Ill CI) . The latter case will naturally lead us to in- 
troduce the concept of electrostatic line tension which 
we will show survives to a large extent in the com- 
plete Poisson-Boltzmann treatment ( Sees . lITTjl and [ITDl) . 
Based on this concept we will then describe the rupture of 
tethered as well as fluid shells and finally construct the 
appropriate stability "phase diagram" of the ruptured 
shells (Sec. IIIII) . We will compare our results with previ- 
ous attempts at elucidation of the stability of spherical 
aggregates and conclude by pointing out the most impor- 
tant consequences of our theory (Sees. ITVl and |V| . 



II. PB LEVEL DESCRIPTION OF 
ELECTROSTATICS OF A PARTIALLY FORMED 
CHARGED SPHERICAL SHELL 

The PB theory of electrostatic interactions in ionic so- 
lutions is based on the mean-field approximation and is 
thus applicable only for sufficiently small charge densi- 
ties on the surface of the macromolecules, low ion va- 
lencies, high medium dielectric constant, or high tem- 
peratures Q. In the case of monovalent salt solutions 
with low surface charge densities these limitations are 
not particularly severe and the results of the PB theory 
are quantitatively correct even when compared to more 
sophisticated approaches [lj| . The limitations of the PB 
approach become practically important only in highly 
charged systems where ion-ion correlation effects begin 
to affect the electrostatic properties of the charged sys- 
tem [ItJ • Though the accuracy of the PB approach can 
be systematically improved by perturbative corrections 
in the ionic fluctuations around the mean-field solution, 
eventually one has to give up the idea of a mean-field de- 
scription altogether taking recourse in a fundamentally 
different reformulation of the electrostatic theory based 
on the concept of strong coupling 3] . We shall not delve 
any further into these reformulations assuming in what 
follows that the solution conditions we are describing are 
far removed from the limit of strong coupling. 



A. Model of a partially formed spherical shell and 
modes of rupture of tethered and fluid shells 

The geometry of the partially formed spherical shell is 
shown in panel a) of Fig. [T] The shape can be charac- 
terized by two parameters: the shell radius (R) and the 
opening angle ($o)- The surface of the shell is assumed 
to be homogeneously charged, so that the surface charge 
density a is a constant. Panels b) and c) of Fig. [T] sketch 
the ruptured shells in the case of tethered and fluid shells, 
respectively. 

We represent the rupture of a charged tethered shell in 
terms of a partition of the intact solid shell into two sepa- 
rate pieces described as two spherical caps with opening 
angles $o and tt — $0, with radii equal to that of the 
complete original shell (R). The surface charge density 
is assumed to remain the same in both pieces of the rup- 
tured shell and equal to the one of the complete shell 

In the case of a fluid shell we assume that the rupture 
starts via formation of a pore that increases its radius as 
the rupture progresses and the pore opens up. During 
the poration the partially opened fluid shell is described 
as a spherical cap of radius R($o)i which is not equal to 
the radius of the original shell before the poration. We 
shall presume that the density of the shell material and 
the shell thickness remain the same during rupture and 
that the material of the fluid shell is conserved during 
poration. This also presumes that the total charge of the 
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shell is conserved, so that the surface charge density does 
not change in the process of poration. 





R 





b) 







FIG. 1: a) The geometry of a ruptured spherical shell, b) The 
mode of rupture of a tethered (solid) shell, c) The mode of 
rupture of a fluid shell. The amount of shell material in both 
modes of rupture is conserved in the process but is partitioned 
between the original shell and its progeny in the case of a 
tethered shell, or is contained completely in the progeny in 
the case of a fluid shell. 



B. Electrostatic potential and free energy of a 
partially formed shell at PB level 

In the following, we shall assume that the bathing so- 
lution of the ruptured shell is an aqueous monovalent 
salt solution (relative permittivity e) whose bulk concen- 
tration equals Co- We describe the electrostatics of the 
partially formed shell in the PB framework [l8|, [l9| . Our 
approach to the problem is the same as that described 
in Ref. l20t however, the geometry is more involved in 
this case. In short, to obtain the electrostatic potential 
(p in the solution that in this case depends on the radial 
coordinate, r, as well as the azimuthal angle, i9, we solve 
the PB equation, 



VV(r, 0) = sinh (pe tp(r, I?) 

ee V 



(1) 



subjected to ^-dependent standard electrostatic bound- 
ary condition on the shell surface of the form 



d<p(r, 



Or 



d(p(r, i3) 



dr 



--R+ 



ee 



for de[d ,TT), 
(2) 

where j3 1 = fcgT with Izb the Boltzmann constant and 
T temperature, e the vacuum permittivity, and e the el- 
ementary charge. We assume an infinitely dilute solution 
of ruptured shells so that there are no interactions either 
between intact shells or between the products of the rup- 
ture. We assume a salt bath in chemical equilibrium with 
the (ruptured) shell and its bathing ionic solution. The 
concentrations of positive and negative ions in the bath 
are equal, (P + = c°_ = Co . 

Once the electrostatic potential is obtained, the free 
energy of the system can be constructed by standard 
methods as detailed in Ref. [2(| The electrostatic part of 



the (Helmholtz) free energy (F e i) can be evaluated in var- 
ious equivalent ways but is most conveniently calculated 
from the charging process 



AS / ip(r = da 



(3) 



where the upper bound of the integration over the surface 
charge density is set by the value of the surface charge 
residing on the surface S [21] . The above form of the 
free energy can be easily simplified in the case of the DH 
linearization of the PB theory by acknowledging the fact 
that the DH electrostatic potential depends linearly on 
the charge density. 



C. Analytical solution of the problem in DH 
approximation 

When the electrostatic potentials in the electrolyte so- 
lution are small, /3eo</? -C 1, the PB equation can be 
linearized leading to the DH equation of the Helmholtz 
type whose solutions are well known. 



1. DH equation, boundary condition, and free energy 

Linearizing the PB equation [Eq. (pj], the DH equa- 
tion is then obtained as 

V 2 <p(r,ti)- K 2 (p(r,d) =0 , (4) 

where k is the inverse DH screening length given by 

= y / ee /2/3e^c . (5) 

The electrostatic problem of a shell with radius R that 
has a hole characterized by opening angle t9q can be posed 
also in slightly different terms - we may treat it as a 
problem of a complete spherical shell but with an az- 
imuthally varying surface charge density, <r(i?). In this 
case the surface charge density assumes a Heaviside-like 
form in terms of the azimuthal angle 



tr(0) = 



, i? < <& 



a , d < -d < n 
The boundary condition in Eq. ([2]) then becomes 



d<p(r, $) 



Or 



dtp(r, i?) 



R=R~ 



Or 



R=R 



+ 



es 



(6) 



(7) 



The potential that is a solution of Eq. (Q| can be 
written as a superposition (22[, 



¥>{r, #) = X] ~7= aiI l+ i(K,r) + biK l+ i(Kr) Pi(cosd) 

1=0 V r " 



(8) 
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where a; and bi are the expansion coefficients, Pi are Leg- 
endre polynomials, and I l+ i and K l+ i are the modified 
spherical Bessel functions of the first and the third kind, 
respectively. 

To obtain the free energy of the shell corresponding 
to the appropriate solution of the DH equation (see Ap- 
pendix [AJ , one needs to evaluate the charging integral 
given by Eq. ([3]). The electrostatic free energy in the 
DH approximation can then be obtained in the form 



/•7T 

F^ H =TvR 2 a / dtf sin <p{R, d) 



yielding finally 



ixR 2 a 2 1 



ee oK 2j-^ (21 + 1)F (1,kR) 



(9) 



(10) 



The functions FaQ, $o) and Ta(l, kR) are defined in Eqs. 
HI) and (1X41). 



Asymptotic expansion and limiting form of free energy 



The numerical evaluation of Eq. (jlOl) does not pose 
significant problems, but the analytical limits that may 
be of interest turn out to be more difficult to come by. 
These were obtained using routines from Mathcmatica 

& 

One is first interested in the limit when $o ~~ ^ 0. This 
represents the complete shell without a hole. The result 



irR 2 <7 2 



lim F^ M = 

»?o->0 EEqK 1 



coth kR 



(11) 



which is the same as Eq. (15) in Ref. [20. 

Some insight can be obtained by analyzing the asymp- 
totic forms of the Bessel functions featuring in the free 
energy- By using the asymptotic expansions, given in 
Refs. [2J and [2a. we obtain 



FT 



nR 2 (7 2 
eeoK 2 



1 oo 

IE 



1=0 



21 + 1 2 



1 



1(1 + 1) 
2p 2 



-O 



1 



,P 
(12) 

where we have introduced p = kR. It turns out that the 
asymptotic expansion is problematic for any finite value 
of p, so that the free energy for an open shell can only 
be obtained by a significantly extended summation of 
the nonasymptotic I series (detailed discussion is given in 
Appendix [Bj) . Since effectively no asymptotic result can 
be derived via a simplistic expansion of the free energy 
in Eq. (ITU)) , one needs to approach the problem by a 
different route. Note here, however, that the zeroth-order 
terms (proportional to p°) of the expansion [Eq. (fl~2j)] 
may well be valid, as they do not depend on the effects 
of the edge and its screening. We shall establish in what 
follows that this is indeed the case. 



Problems with the range of validity of the asymptotic 
form of the electrostatic free energy obtained directly 
from the solution of the DH equation lead us to an al- 
ternative formulation that would avoid the pitfalls of the 
infinite summation due to the sharp edges of the partially 
opened shell. It is based on the simple equivalence be- 
tween two ways of calculating the electrostatic energy of 
the system: it can be thought of as a result of the charg- 
ing of the system or as a result of the direct interactions 
in the system. Both viewpoints are equivalent plj . 

By taking now the second approach and summing the 
interactions over the surface of the incomplete shell (see 
Eq. (|C2[) and the derivation in the Appendix [C]) , one 
obtains the form of the free energy in the limit when 
kR 3> 1 but remains finite: 



lim FS H 



7rR 2 o~ 2 

eeoK 



1 + cos do f 



kR 



sin -dg 



(13) 



where / is a numerical constant, / ~ 0.12. The second 
term in Eq. (fT5)) arises from the fact that the charges 
on the shell close to the open shell edge lack some of 
the neighboring charges with respect to the charges away 
from the edge. The parameter that determines the effec- 
tive distance of a particular charge from the shell edge is 
the screening length, 1/k. Since the length of the open 
boundary of the shell is 



C = 27ri?sini? , 



(14) 



the second term is obviously proportional to C and can be 
interpreted as the interaction renormalized line tension 
energy of the edge. The interaction renormalization of 
mesoscopic material properties, line tension in this case, 
is a ubiquitous feature of molecular interactions in matter 
[2(| and its appearance should also not come as too big 
a surprise in the present context. 

Note, however, that since the result was derived in the 
limit when kR 3> 1 but still finite, the second term is 
significantly smaller from the first one as the contribution 
of the edge to the total electrostatic energy of the shell is 
small in this limit. For e.g. viruses, the condition kR ^> 1 
is typically well obeyed in physiological conditions 0, 

Now that we have found an appropriate form of the 
electrostatic free energy in the limit 00 > kR > 1 we can 
check how well the exact (non-asymptotic) DH solution 
[Eq. [TUJ fits to it. One should recall that we were not 
able to obtain this form of the free energy by merely 
analyzing the asymptotic form of the DH free energy. In 
accordance with Eq. (IT51) . we take the fitting form 

f>«(0o; a,b) = a( ^ + C ° B ^ - frsinflo^ (15) 

where a and b are numerical constants to be determined 
by the fit to Eq. flU}]). The fit is shown in Fig. H We 
find a non-zero value for b of the form 



kR 



where 6n ss 0.161 



(16) 



5 



This is completely consistent with the approximate form 
in Eq. (fT3")l , and the value of / that we obtained from 
the approximate analysis in kR ^> 1 regime is not far off 
from the value of bo- 




FIG. 2: The results of the fit of Eq. (JT5J) to the DH solution 
(pluses). The full line shows the bo/p behavior with bo ~ 
0.161 . The empty symbols show the behavior of the fit to 
Eq. ([151) of the full PB solution calculated for R = 10.04 
nm, for different surface charge densities: a = 0.1 eo/nm 2 
(squares), 0.4 eo/nm 2 (circles), 0.8 eo/nm 2 (triangles), and 
1.2 eo/nm 2 (diamonds). Dashed lines are guides to the eye. 
The numerical solution of the PB solution become less reliable 
with increasing p as the errorbars indicate. A noticeable effect 
in PB case is the dependence of the line tension coefficient b 
on the surface charge density. 

The deviations from the 1/p behavior are apparent in 
the regime where p < 1 , where in fact we do not ob- 
serve the part with sini?o dependence predicted by Eq. 
(|15[) and the line tension renormalization concept is of 
limited use there. In this regime, the screening from the 
inside of the shell is incomplete since the volume of the 
shell interior is too small to contain all the salt ions re- 
quired to screen the shell charge, i.e. R is too small since 
p ~ 1 . Thus, the line tension renormalization is differ- 
ent from the form implied by Eq. (fT3")) . This effect only 
becomes important for small angles i?o when there is a 
clear separation between interior and exterior. For larger 
values of t?o , when the shell looks more like a disk than a 
sphere, the separation of space into interior and exterior 
does not make sense and sint?o behavior holds for large 
enough angles, even when p ~ 1 . 



D. Solution of nonlinear PB equation 

We now solve the full nonlinear PB equation [Eq. (JlJ] 
and evaluate the corresponding electrostatic free energy 
[Eq. ©]. 

Due to the axial symmetry, the problem is effectively 
two-dimensional, and the solution can be sought in the 



domain (r, <E [0,oo) x [0,7r). To solve Eq. (P), we 
first discretize the domain and the PB equation for the 
electrostatic potential using finite differences. Typically, 
we used about 1200 points in the radial and 500 in the 
polar direction in the nonlinear regime of parameters. 
The largest radial point has to be chosen such that the 
potential at this point is very nearly close to zero. The 
band matrix system obtained in this way is solved using 
a damped nonlinear Newton method [28(, iterated until 
desired convergence is reached. 

Figure [3] shows the numerical solution of Eq. ([1} for 
the dimcnsionlcss electrostatic potential f3eo<f- 






x[nm] 

FIG. 3: The solution of the PB equation for the dimensionless 
electrostatic potential /3eo<fi . The shell radius is R = 10.04 
nm, and the opening angle $o = 0.2357T . The shell surface 
charge density is a = 0.4 eo/nm 2 , and the bulk concentration 
of salt ions is co = 10 mM. 

Figure^] shows the dependence of the electrostatic free 
energy of the system on the opening angle of the shell. In 
the limit when i3q ^ , the problem reduces to the one 
of a closed shell studied in Ref. [2(| these results are also 
shown in Fig. 01 The observed agreement nicely parallels 
the one obtained in the DH limit as well, since we have 
previously demonstrated that in the limit when z?o 
the analytical DH result [Eq. (fTU)) ] reduces to the form 
previously derived in Ref. [2fJ. We also compare the DH 
results obtained from Eq. (fTU|) and shown in Fig. 2] with 
the numerically exact PB results. 

We now examine the concept of electrostatic renormal- 
ization of the line tension that we introduced before at 
the DH level. The question is whether and how the con- 
cept survives in the numerically exact solution of the PB 
theory in the region of parameters where the DH approx- 
imation does not apply. To test the usefulness of the line 
tension concept even outside the DH limit, we have com- 
pared the angular dependence of the electrostatic free en- 
ergy obtained from a numerical solution of the complete 
PB equation to the form of Eq. (|15l) . The result of a 
numerical solution of the PB equation is shown in Fig. [5j 
and the values of line tension coefficients b obtained from 
fits to Eq. (|T5]) are summarized in Fig. [5] 

It is obvious from Fig. [5] that the line tension depen- 
dence describes the behavior of the electrostatic free en- 
ergy quite well even in the case of complete nonlinear 
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FIG. 4: Electrostatic free energy of the system as a function 
of the opening angle, for a capsid of radius R — 10.04 nm 
and several values of surface charge densities and salt concen- 
trations, as denoted in the body of the figure. The dashed 
lines show the solution obtained in the DH approximation; 
the black dashed line aligns almost completely with the nu- 
merical PB solution. The star symbols show the PB results 
for a closed capsid (f?o = 0) taken from Ref. [20l . 

PB theory. The numerical data, even in the case of very 
small salt concentrations (e.g. 1 mM) where the DH ap- 
proximation does not apply, can be fitted to the general 
form of Eq. (IT5|) . although the electrostatic contribution 
to the line tension coefficient, b, may be quite different 
from the DH prediction. As expected, the numerical data 
approach the DH prediction as the salt concentration in- 
creases, as is usual since the electrostatic potentials in the 
solution decrease [2(| ■ Contrary to the prediction of the 
DH theory [Eq. l[T5|) ] the line tension coefficient shows 
dependence on the surface charge density (see Fig. [2]), 
particularly in regimes where the DH limit docs not hold. 

III. RUPTURE OF SHELLS 

Regardless of the smallness of the electrostatic con- 
tribution to the edge tension energy, particularly when 
p^$> 1 , it is precisely this term that induces the shell rup- 
ture. One notices that electrostatics favors the opening of 
the edge. It is the opposing attractive interactions of the 
vdW/hydrophobic type that keep the shell together. We 
shall assume that the attractive interactions are of short 
range, so that the attractive part of the cohesive energy 




0.0 0.2 0.4 0.6 0.8 1.0 
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FIG. 5: Fit of the angular dependence of the solution of the 
PB theory to the form Ffui'&o) = a[(l + cosi?o)/2 — ftsintfo], 
for a shell of radius R = 10.04 nm and surface charge den- 
sity a — 0.4 eo/nm 2 . Salt concentration is 10 mM. The 
inset shows a separate fit of — fesintfo to the numerical so- 
lution, from which the cosine dependence was removed, i.e. 
F n um.("&o) / F nU m{0) — (1 + cos$o)/2 . Note that in the latter 
case, a slight effect of numerical error can be discerned. 



can be approximately associated with "bonds" between 
the molecular constituents of the shell. Hydrophobic and 
van der Waals interactions could be conceived to be of 
such nature. These soft molecular bonds are broken upon 
rupture, so that the part of the cohesive energy that is 
lost in the process is proportional to the edge length. 

This cohesive energy change is different for the case of 
tethered as opposed to fluid shells. In the former case, 
two partially opened shells are created upon rupture, 
whereas in the latter case the rupture creates a single 
partially opened shell. Also in the case of the fluid shell 
rupture one has to properly include the change in the 
curvature energy upon rupture, since the pore changes 
the radius of curvature of the shell (see Fig. [TJ . We thus 
treat the cases of tethered and fluid shells separately. 



A. Rupture of tethered shells 

Tethered shells with a spherical geometry show a more 
or less pronounced faceting, but in what follows we will 
confine ourselves to the ideal sphere, ignoring the effects 
of (icosahedral) faceting [2^, [3(| for the sake of calcula- 
tional feasibility. 
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We denote the line tension parameter arising solely 
from the attractive interactions by 7, so that the total 
energy required to form an edge of length C (or to sepa- 
rate a complete shell into two spherical caps with opening 
angles $0 and n — # ) is 



Fedge 



7' 



boa 2 



C . 



(17) 



Note here that F ec i ge contains the electrostatic contribu- 
tion to the edge energy of the incomplete shell twice, 
because upon rupture, two incomplete shells are formed, 
both having an open edge (note also that, since both 
spherical caps have the same curvature, we do not need 
to include the change in the bending energy here). 

When the quantity in the brackets (the edge tension 
parameter renormalized by electrostatic interaction) in 
the rhs of Eq. (jTTJ) becomes zero, the rupture sponta- 
neously takes place. Since the energy gained is propor- 
tional to the edge length, one can expect the breakage 
of all the bonds, a complete disintegration of the teth- 
ered shell into its (macro)molecular constituents. In the 
context of our model, a particular macromolecule may 
be imagined as a small spherical cap with a certain edge 
length. The rupture of a sphere would thus correspond 
to disintegration of all the small spherical caps from the 
assembled structure. Such final state would have the low- 
est free energy of all others that could be formed by as- 
sembling the small spherical caps (macromolecules) in 
several larger cap-like aggregates (ruptured parts of the 
shell); this conclusion holds as long the concept of line 
tension makes sense. This finding enables us to draw a 
diagram of the tethered shell stability, i.e. to draw a line 
separating the stable shells from the ruptured ones in the 
(j — Co plane (with 7 assumed to be constant and inde- 
pendent of <7 and Co). This phase diagram is shown in 
Fig. E 

There arises a question on the reliability of the insta- 
bility line obtained in this way. Namely, the DH approx- 
imation is reliable only when the electrostatic potential 
if in the solution is small, i.e. when fie$Lp <C 1. We 
have shown, however, that the concept of the electrostatic 
renormalization of the line tension carries over from the 
DH limit onto the full PB theory almost unchanged. The 
most important difference between the DH and the com- 
plete PB theory is that the electrostatic renormalization 
of the line tension is an effect significantly smaller than 
predicted by the DH approximation (see Fig. [5]). This 
is also the reason that in the (co,cr) phase diagram the 
shells are stable in a larger regime than predicted by the 
DH limit (see Fig. [BJ. 



B. Poration (rupture) of fluid membrane shells 

As already discussed, fluid membrane shell can rupture 
without loss of material by simply rearranging the con- 
stituent (lipid) molecules so that the hole opens up. It 
is known that the vesicles may form holes in spite of the 
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FIG. 6: Symbols: Instability line of tethered shells towards 
rupture obtained by numerical PB calculation with 7 = 1 
fcsT/nm and shell radius R =10.04 nm. Dashed line: Debye- 
Hiickel approximation, a = eo\/2(3jco/bo- In PB case, the 
points (a, Co) of zero edge energy were obtained by interpola- 
tion of numerical results for several surface charge densities a 
and salt concentrations Co, hence the errorbars. 



line tension (exposure of hydrophobic parts of the con- 
stituent molecules) if their surface is charged |3lj] . Ex- 
perimentally, such pores have been observed for instance 
in the so-called red blood cell ghosts at low salt concen- 
trations [32| . 

Here we briefly investigate the possibility that the 
charge on the vesicle induces a different kind of surface 
deformation. Instead of vesicle poration, one may think 
of a situation where the vesicle separates into two smaller, 
completely closed vesicles, so that the edge does not form 
in the process. This scenario is similar to the Rayleigh 
instability of charged droplets Q discussed in the Intro- 
duction. However, since the electrostatic energy of the 
charged shell in the DH regime (and when kR 3> 1) is 
proportional to 20] a 2 R 2 , and in the process of transfor- 
mation into two smaller closed vesicles the area remains 
the same, there is no gain in the electrostatic energy in 
the final state of the system. This type of instability is 
of importance only in the poor screening, near-Coulomb 
regime, where the electrostatic energy of the vesicle is 
proportional to the third power of the vesicle radius 20] . 
We shall neglect these effects in the following and con- 
centrate exclusively on the charge induced poration. 

It can be envisaged that the vesicle with a hole may 
have quite complex geometry (33|, but we shall approxi- 
mate it by a sphere with a missing cap as in the previous 
sections. Due to the way the pore forms we consider 
the sequence of shapes that have the same number of 
molecules and thus the same area, but different radii, or 
equivalently, the opening angles of the pore. Such a se- 
quence would thus start from a fully closed vesicle and 
open to a flattened-out disk. We parametrize this se- 
quence of shapes with an opening angle i?o of the formed 
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pore. 

For the total area of the vesicle with a pore to equal 
the area of a vesicle without pore the radii of the shapes 
in the sequence must grow with the opening of the pore, 
R = i?($o) . Given a closed vesicle with radius Rq and 
area A = AttRq , some elementary geometry then gives 



iZ(0o) = Ro 



1 + COS $i 







(18) 



and similarly for the p($o) — K R{$o) > where we now 
denote po — kRq . 

We need to consider the influence of the elastic bend- 
ing energy of the vesicle as well, since the curvature of 
the shell changes as the pore opens. In the classical cur- 
vature model, the local elastic bending energy is in the 
lowest order composed of two terms containing mean and 
Gaussian curvature [HJ . With the two curvatures we as- 
sociate two elastic constants, n and kq, called bending 
rigidity and Gaussian bending rigidity, respectively. As 
we are interested only in a piecewise spherical surface, 
both radii of curvature are the same. Using i?($o) de- 
fined above we see that the curvature free energy can be 
written as 1351 



Fbend — 1"kK ■ 



1 + COS 1? 







(19) 



Here the Gaussian and the bending elastic moduli have 
been combined into an effective bending rigidity param- 
eter of the vesicle, JC = 4(k + kg/ 2) ■ In the equation 
above we have assumed that the spontaneous curvature 
of the shell is zero, although the case with non- vanishing 
value of the spontaneous curvature can be easily studied 
along the lines that we present in the following. 

We now derive the complete free energy F to tai = 
jC+Fbend+F® 11 including thus the non-electrostatic line 
tension energy, the bending energy, and the electrostatic 
energy. The latter we derive first within the DH approx- 
imation confining ourselves to the limit kR^q) 3> 1. In 
this case the electrostatic part of the free energy can be 
derived in the DH limit by inserting R(i)o) into Eq. (p~3|) 
and obtaining for the total free energy 

F tota i = ^^+F total {d ) , (20) 



E£qK 



where 



Ftotai (#o) 1+COS^o 



2tt/C 



V F ( 1 - iSL ) V2(l-costf ) ■ (21) 



Here we introduced the vesiculation index Vf and the 
electrostatic line tension 7 e ; as 

VF = — =2^C and 7ei = W' (22) 



with Cq = 2-kRq . Note that Vp is defined as half the value 
introduced in Ref. l36l . The only parameters determining 
relative contributions of different free energies are thus 
Vf , the ratio "f e i/l , and po ■ 

The specific rewrite of the energetics of the shell in 
the DH limit with p » 1 reduces the problem to the 
one studied by Fromherz in Ref. 36], where only the line 
tension and curvature contributions were considered. In 
our case, however, we have an additional possibility of 
the line tension energy being both positive and negative 
as a consequence of the electrostatic effects. Without 
electrostatic interactions it is always costly to form edges 
due to residual attractive interactions giving rise to non- 
electrostatic line tension energy, 7. With electrostatic 
interactions included, it may be favorable to form holes. 
This happens when (1 — "f e i/"f) < 0. For pores in the 
vesicle to open up spontaneously, the free energy has to 
fall off as a function of the opening angle. 

It is convenient to introduce the ratio of energy param- 
eters, 



p= V F (l -Jei/l) ■ 



(23) 



The total shell free energy can be written as a functional 
whose universal behavior depends on the parameters p 
and t?o • In Fig. [3 we plot the function 



1 dF 



total 



sin?9n 



P 



2vr/C <9tfn 



sinf 



(24) 



which allows us to separate regions where it is larger 
or smaller than zero (note that JC > 0). The "phase 
diagram" shows the type of information already studied 
by Fromherz, but modified by electrostatic effects in our 
case. One can see three different regimes: The first one 
takes place for when p < 0, where even pores of the 
smallest size are unstable and they thus spontaneously 
form. For < p < 1 , only the pores of certain critical 
radius will spontaneously grow, otherwise they will close. 
The critical pore opening angle (i.e. the critical pore 
radius) for given p is at the point where the white region 
meets the gray regions. For p > 1, there are no stable 
pores. The second regime, where < p < 1, can be 
divided in two parts: in the first one (0 < p < 0.5) the 
disk state has lower free energy than the sphere state, 
and in the second one (0.5 < p < 1) it is the other way 
around; both regimes are indicated in Fig. [7J 

The PB approach only shifts the borders in the DH 
phase diagram, which remains qualitatively the same, i.e. 
we observe no additional effects, metastable pores in par- 
ticular. This is a consequence of the validity of the line 
tension concept even in the PB limit. As we have shown 
(see Fig. [5]), the PB solution only predicts a lower value 
of line tension, the more so the larger the surface charge 
density. This effect cannot change the salient features of 
the pore energetics discussed previously at the DH level. 

We note in passing that the formation of pores could 
be also described as pore nucleation, emphasizing in this 
way the similarity of the process to the physics of the 
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FIG. 7: The "phase diagram" of the ruptured capsids that 
shows the sign of dF tota il d"&v> (gray regions). In light-gray 
(dark-gray) regions, the fully opened, disk state of the mem- 
brane has lower (higher) free energy than the spherical state. 
In the white regions, the free energy always increases with $o 
and the vesicle is stable towards pore formation. 



classical theory of nucleation [37J . As is the case there, 
we also have a critical size of the pore required in order 
to start the pore nucleation process, which extends the 
pore to the maximal possible size, leaving the fluid shell 
in the disk state. 

Let us now relate the results obtained for fluid shells 
to those previously obtained for the tethered shells. In 
the case of tethered shells, the energy contribution of 
bending rigidity of the shell was of no importance for 
the rupture mode. If we put K, = 0, then the phase 
diagram shown in Fig. [7] separates in two regions only, 
the one where 1 — "f e ih > 0, p — > +oo, and the shell 
is stable towards pore formation, and the other where 
1 ~ 7e//7 < 0, p — > — oo, and pores spontaneously form. 
Without consideration of the bending rigidity, the region 
of parameters where the pores of critical radius grow and 
smaller ones close does not exist. The stability line is 
given by the condition 1 — -f e i /-f = which gives 



bo * 2 _ 

7 - -w 2 = 7 - 7ei = . 

2 eenn 



This is the same critical relation as the one following 
from Eq. (IT71) up to the multiplicative constant of 2 
which arises due to two electrostatically renormalized 
open edges that exist in the case of tethered shells. 
Therefore, when K. — 0, the discussion regarding the sta- 
bility of fluid shells towards poration is quite similar as 
in the case of tethered shells discussed previously. As in 
that case, in a region of charge densities and salt concen- 
trations above the full line in Fig. [5] (er oc ^/cq] large cr's, 
small cq's), the shells will spontaneously rupture. 



IV. APPLICATION OF RESULTS TO VIRUSES 
AND VESICLES 

An infinitely thin, spherical, and uniformly charged 
tethered shell may seem like an oversimplified model of 
an icosahedral virus capsid. Indeed, icosahedral viruses 
are not perfectly spherical [2!| |3(| , their capsids are sev- 
eral nanometers thick, and the charges on their con- 
stituent proteins are not uniformly distributed on the 
capsid surface but show variations consistent with their 
icosa(delta)hedral symmetry [38|. One could formulate a 
more complicated electrostatic model of a virus capsid, 
as has been done in the literature [2(| [H| . One should 
remember, however, that our primary interest is in the 
capsid rupture and the energetics of edge formation. The 
almost analytical type of insight in this matter that we 
detailed in the previous sections would be difficult to ob- 
tain with more complex capsid models that would need 
to rely completely on numerical methods [38j . 

The assumption of vanishing thickness of the shell may 
appear as the most worrying feature of our model when 
applied to real viruses. Indeed, typical capsid thickness 
and the Debye-Hiickel screening length at physiological 
ionic strength are quite comparable. As the charge den- 
sity in general varies across the capsid thickness [39| , one 
may infer that this could invalidate the application of our 
model to real viruses. This is not the case since the pro- 
tein material of the capsid cannot be treated as totally 
permeable to salt ions, so that there is no screening of 
charge across the capsid thickness. Instead, the protein 
material should be treated as a medium with small rel- 
ative dielectric constant, as has been done e.g. in Refs. 
l2fj| and H3- This has been studied in some detail pre- 
viously in Refill, where it was demonstrated that the 
most important features of the spherical shell energetics 
are already contained in the model of an infinitely thin 
charged shell. This also includes the analytical behavior 
of the electrostatic energy, i.e. its scaling with the shell 
radius and the surface charge density. 

Our calculations predict an instability of a charged 
tethered shell that depends on its surface charge density 
and the concentration of monovalent salt in the solution 
into which it is immersed. The tethered shells become 
unstable cither with (i) increase of the surface charge 
density, or (ii) with decrease of the salt concentration. 
This holds both in the DH and PB approaches to the 
problem as can be seen in Fig. [6] These results should 
apply to viruses, but to the best of our knowledge, there 
are no systematic studies of empty capsid stability with 
respect to changes in pH and salinity. There are studies 
of the conditions for the self-assembly of empty capsids 
[ill - lii} . These are not directly applicable to the problem 
of interest to us, since the assembly necessarily involves 
change in the entropy of dissolved proteins (monomers) 
(35l l4a | . Assuming that the capsid disassembly and as- 
sembly are two faces of the same thermodynamical equi- 
librium process, one finds that the capsid stability de- 
pends on the critical concentration of proteins which is 
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related to the free energy gain (loss) upon assembly (dis- 
assembly) [H, HBj]. Thus, in sufficiently low concentra- 
tions, the capsids should be unstable, regardless of the 
strength of the protein contacts in the capsid. However, 
it has been demonstrated that the disassembly of hepati- 
tis B virus capsids proceeds exceedingly slowly and that 
the capsids are stable for days in the low-concentration 
conditions where the disassembl y sh ould take place [46| . 
This has led the authors of Ref. |46| to conclude that the 
capsid disassembly and assembly reactions are not in the 
equilibrium. This is why our approach should be of use, 
as it concentrates on the enthalpic contribution to the 
shell energetics, and the disassembly induced in such a 
way should depend on the changes in ionic strength and 
pH of the solution and not on the concentration of as- 
sembled capsids. The long-term stability of capsids in 
low-concentration conditions in fact enables experiments 
that could check our predictions. 

The concept of line tension energy has been shown to 
play an important role in the nucleation models of capsid 
assembly [H, [45| . While it does not importantly influ- 
ence the distribution of capsids and protein clusters in 
the assembled, equilibrium state [35[, its presence does 
lead to the appearance of a barrier in the assembly free 
energy, resulting in sigmoidal kinetics of capsid assembly, 
i.e. a characteristic lag time that depends on the barrier 
height 45] . We have demonstrated that the line tension 
is necessarily modified by electrostatic effects, and thus 
the lag time should vary with the salinity and the surface 
charge density, which is an effect that has only partially 
been accounted for in the existing capsid nucleation the- 
ories [45|, |47|. Our study also demonstrates that some 
care must be exercised when applying the approximate 
DH results, since these may be unreliable in the regime of 
surface charge densities and salinities relevant to viruses 
(see Fig. [5]) , in accordance with the findings in Ref. H^. 

In the case of liquid shells, our model is best applied to 
rupture of charged and sufficiently thin vesicles - in that 
case the surface charge we use corresponds approximately 
to the monopole surface charge of the vesicle. The model 
is also useful, however, when the vesicle thickness cannot 
be neglected on the scale of the screening length, as has 
already been discussed in the case of viruses (see also 
Refs. [48|that represent the spherical vesicle as two shells, 
bilayer of charge). In the simpler variant of the fluid 
shell model considered by Fromherz, no stable minima 
corresponding to intermediate shapes occur 36]. The 
same is seen from our model, in spite of the fact that 
in addition we also consider the electrostatic effects - 
the only shapes that minimize the free energy are the 
closed vesicle and the flattened disk. This would lead 
directly to an all-or-none transition in the assembly of 
the shells. Indeed it has been pointed out at various 
times that e.g. the icosahedral capsid assembly can be 
construed in analogy to a continuous phase transition 
[HI, , thus requiring no latent heat in the process of 
assembly that furthermore proceeds without a nucleation 
process. 



While Fromherz studied a model of the vesicle that 
contained the non-electrostatic line tension energy and 
bending energy contributions, Betterton and Brenner [3lJ 
studied the vesicle poration neglecting the bending en- 
ergy contribution, but accounting (at the DH level) for 
the electrostatic and non-electrostatic contributions to 
the line tension energy. Our approach is more general 
than the two previous approaches, and contains both of 
them as limiting cases. 

In the study of Betterton and Brenner, the vesicle was 
modeled as a flat membrane of constant area and charge 
and the electrostatics was considered in the DH limit 
only. For this approximation to hold, the pore radius 
should be small compared to the curvature of the vesicle. 
It was shown that in such a model there exists a narrow 
regime where metastable pores occur [3l|. We did not 
observe these effects neither in the DH approximation 
nor in the full PB solution, even when kR < 1, i.e. in the 
regime where the line tension renormalization concept is 
of limited use. To stabilize the pore in fluid membrane, 
mechanisms other than those considered here are needed. 
These may include adsorption of multivalent ions [50] or 



dynamical stabilization of pores [51 



Although our primary interest is in stability of charged 
shells, our study offers interesting insight in their possible 
modes of growth. If we assume that the spontaneous cur- 
vature of the shell is zero (as we did), one may model the 
shell growth by the addition of subunits to the initially 
flat patch. Our model predicts that a sufficiently large 
planar piece of shell material will spontaneously close into 
a sphere when p > 1, see Eq. ([23]) (the shells of smaller 
radius may form but these require thermal excitations in 
order to cross a free energy barrier, see Fig. [7])- Again we 
should emphasize that no intermediate shapes between a 
flat patch and a complete spherical shell exist. The ra- 
dius of thus formed spherical shell is 



Rv > 



K 



7-Jel 



(26) 



This inequality shows that the charged shell radius in 
general depends on its surface charge density and the 
salinity of the solution. We apply Eq. (|26l) to viruses 
using 7 ~ 2 - 8 k B T [H, K ~ 60 k B T [H, c = 100 
mM, and o ~ 0.5 eo/nm 2 , which yields Rv ~ 10 — 30 ran. 
This is indeed a range of radii typical for viruses. We note 
here that the electrostatic renormalization of line tension 
is quite small, only ~ 0.2 k B T, and has thus a small ef- 
fect on the radius of the virus. The electrostatic effects 
become more important in solutions of smaller salinities 
which in principle enables formation of larger viruses. 
This mode of growth is different from nucleation models 
of virus capsids that assume a fixed curvature of the grow- 
ing shell [H, HH , but bears some conceptual similarity to 
other proposed modes of capsid growth [l3| . As in pre- 
vious nucleation models, it is still difficult to explain the 
characteristic Caspar-Klug type of order of proteins in a 
complete capsid that is the hallmark of simple viruses. 
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SUMMARY AND CONCLUSION 



Appendix A: Solution of DH equation 



We have performed a detailed study of the contribu- 
tion of electrostatic interactions to the overall energetics 
of an open charged spherical shell immersed in a solution 
of monovalent salt. The results we obtained can be used 
to assess the stability of charged vesicles and viruses. In 
addition, our study also offers some fresh insight regard- 
ing the growth of charged shells and predicts a lower limit 
on the charged shell radius that depends on the physi- 
cal properties of the shell material as demonstrated in 
Eq. (1261) . The alternative description we propose pro- 
vides a characteristic radius of the charged shell without 
invoking the spontaneous curvature of the shell material 

We have proven by a detailed analysis of the numerical 
solutions of the PB equation as well as by an analytical 
solution of the DH approximation that electrostatic ef- 
fects lead to only two stable shapes of an open spherical 
shell, without any stable intermediates, that minimize 
the total free energy. These shapes are the closed vesicle 
and the flattened disk. 

Concerning the rupture of vesicles, our study improves 
on previous ones in several important aspects. When 
considering electrostatic effects we account for the full 
spherical geometry and shape morphology of the vesi- 
cle. Furthermore, for such a geometry (including a pore) 
we have been able to analytically solve the problem at 
the DH level of approximation. This should be of use 
in a broader context of the electrostatic self-energy of 
complicated shapes [see e.g Ref. [l| and references therein] 
that do not lend themselves to an easy explicit solution 
of the DH equation but can be analyzed approximately 
by the methods introduced in this work. Finally, the 
numerical procedure that we developed enabled us to 
study the nonlinear effects arising in the full Poisson- 
Boltzmann approach to the problem by showing that the 
insight given by the DH solution persists also in the non- 
linear PB regime. By combining all of this, we identified 
possible pitfalls in the previous studies that have almost 
exclusively treated the similar class of problems at the 
DH level [27L ED, El. 



In order to solve Eq. [4] in the required geometry, we 
shall partition the space into the internal and external 
region of the shell, defined by r < R (region I) and r > 
R (region II), respectively. Since K l+ i(nr) —> oo when 
r — > and I l+ i(nr) — > oo when r — > oo, the solutions 
inside and outside the shell are given by 



tpi(r,0) = 5^j4,*i(«»Oi}(cofl0) , 

1=0 
oo 

<Pii{r,6) = J2 B i k i( Kr ) p i( cos #) > ( A1 ) 



1=0 



and h = 



where we have introduced ii = 
K l+ i/^¥ (seee.g. Ref. S3). 

Applying the boundary condition in Eq. ([7]) and re- 
quiring also that the potential be continuous at r = R, 
we get two equations that relate the expansion coeffi- 
cients Ai and Bi which enables one to completely solve 
the problem. 

By using both boundary conditions together with el- 
ementary properties of modified spherical Bessel func- 
tions 12211 . one can derive 



_o_ 1 TsMo) 

2eeoK ii(nR) Fo(1,kR) 



and B 



ii{nR) 
HkR)' 



where 



and 



Tq (I, kR) = 



_ Ii+3/2(kR) Ki +3/2 (kR) 



Ii +1/2 (kR) K 1+1/2 {kR) 



{•COS $q 

(21 + 1) I Pi{x)dx 



At , 
(A2) 

(A3) 



P l+1 (cos tf ) - Pi-i (cos ■do) . (A4) 



While deriving the above identity we have used well 
known properties of Legendre polynomials (22|. This 
completes the solution and enables one to write the po- 
tential in the whole space, and in particular the potential 
on the shell surface itself, ip(r = R,$). 
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From Eq. (|12[) we can deduce the expansion coefficients 
for the free energy in the limit p — > oo , 



lim FS H = 



ttR 2 (j 2 1 + cos$ 
eson 2 



(Bl) 



Due to the way the modified Bessel functions of the first 
and third kind appear in the summation, the first order 
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correction in 1/p cancels out and no 1/p term is obtained 
in this limit, but there may be a term in 1/p 2 . The 
former conclusion is consistent with the general analysis 
of Duplantier [53| but the latter is not. 

In the limit of a closed sphere, do — > 0, the asymptotic 
expansion can be summed exactly, yielding 



lim F 



ttR 2 <t 2 

EEqK 



1 + 



(B2) 



[compare this with Eq. (fTTj)]. From this we see that 
indeed there are no first or second order contributions for 
the case of a closed spherical shell in the limit p — > oo , in 
perfect accordance with Duplantier |53j ] . 

Since the only term in the free energy with explicit 
p dependence is the one containing the ratios of Bcsscl 
functions kR)], the validity of the asymptotic ex- 

pansion is determined by the order of the Bessel functions 
I . An estimate of the range of validity of the expansion 
is given in Ref. [25| in the form 



1 . 



(B3) 



Thus, for I much larger than one, the region of validity 
of the asymptotic expansion can be very far out (p 3> 50 
for / = 10 , for instance). 

In addition, the sum over / in the expression for the 
free energy also includes the Z-dependent factors contin- 
gent on Legendre polynomials of orders I + 1 and I — 1 as 
well as the argument cos "do , which stems from the angu- 
lar dependence due to the finite opening of the spherical 
shell. In principle this part of the free energy still has 
to be summed to infinity, regardless of where in the p 
regime we are. 

Let us now try putting all this together in a final as- 
sessment of the validity of the asymptotic expansion of 
the free energy. Since we are dealing with a DH type of 
problem, the sharpness of the edge of the partially opened 
shell is defined on the scale of the screening length 1/k. 
If the sharpness of the edge is on the scale bigger than 
the screening length, the screening will smoothen out all 
the effects of the edge. In analytical terms the azimuthal 
sharpness of the edge must be given by 



R 



1 

P 



(B4) 



in order not to be washed away by screening. How many 
Legendre polynomials (I's) must one sum up in order to 
get this azimuthal sharpness? About 



Aw 



7T/9 



(B5) 



Thus as p — > oo , so does the maximal number of sum- 
mation terms in the free energy, l max — > oo. This means 
that when considering the validity of the asymptotic ex- 
pansion, the lmax is additionally constrained by the fact 
that it needs to represent the edge correctly, so that one 
gets 



8p > 4(irp) 2 



1 



(B6) 



This is the condition for the applicability of the asymp- 
totic expansion to the problem of a partially opened 
spherical shell that has a hole with a sharp edge. It 
is clear that this condition can not be satisfied for any 
value of p in the asymptotic regime (p> 1). 



Appendix C: DH free energy obtained by summing 
interactions between screened charges 

Within the range of applicability of the DH theory the 
free energy of fixed external charges, such as those resid- 
ing on the incomplete spherical shell, can be calculated 
via a direct summation of the screened interaction free 
energies [2l|. If the fixed external charge is described 
with the charge density po (r) then the interaction energy 
of the charges is given by 



F. 



DH 



d 3 rdV po(r)G(r-r')po(r') 



(CI) 



where Q(r— r') is the Green's function of the DH equation 
in empty space. The equivalence of the forms of Eqs. 
(|Cip and © is exact as long as there are no dielectric 
boundaries in the system. If there are, the form of the 
Green's function is modified but the equivalence of the 
two forms remains valid. 

For point-like charges with po{r) = J^i qi5 3 (r — ri) this 
leads to the following form of the electrostatic free energy 
obtained simply by summing the screened electrostatic 
interaction of these charges: 



r „] — 



1 1 



2 4tT££ 



qiq 3 exp( 

i,3 



-KVij) 



(C2) 



where ru = |r< — r,-|. A marked advantage of the repre- 



sentation in Eq. (|C2[) with respect to the form of Eq. (fT0j) 
is that the effective interaction potential or equivalently 
the Green's function entering the former is short-ranged, 
i.e. cut-off on the length scale of 1/k. 

To obtain the electrostatic free energy of the incom- 
plete shell in the DH regime, it is convenient to imagine 
the shell to be populated by the set of point charges of 
appropriate surface density, so that the surface charge 
density is a , and apply the sum in Eq. (|C2I) . Since the 
summation is effectively cut on the scale of , there are 
two different sets of effective point charges on the capsid 
- ones that are close to the shell edge, so that their dis- 
tance from the edge is smaller than k _1 (region denoted 
by 1 in Fig. |8|) , and the others that are located farther 
away from the edge (region denoted by 2 in Fig. [5]). 

The charges in region 2 are surrounded by neighbors 
from all sides, i.e. they are not under the influence of 
the edge. The charges in region 1, on the other hand, 
are influenced by the presence of the edge, as it is lo- 
cated within the salt screening distance. Evaluation of 
Eq. (|C2p may seem to still pose a problem, since we have 
to sum up the interactions on the curved surface of the 
sphere. These summations can be converted to integrals, 



13 




FIG. 8: Different regions in an incomplete shells, selected on 
the criterion of their distance from the edge. An isolated 
point charge is shown in region 2 with the surrounding that 
contributes to its electrostatic energy. 



if we note that in the limit when the screening distance 
is quite small, nR 3> 1 , the important surrounding of 
a particular point charge is practically fiat - this is due 
to the fact that the spherical cap is practically a piece 
of plane when its opening angle is very small. In our 
case the opening angle of a spherical cap surrounding a 
particular point charge is (rei?) -1 , and the assumption of 
planarity is thus satisfied in the limit we are interested 
in, and that, in particular, applies well to viruses [20l[27|. 

For the electrostatic energy of a "patch" of infinitesi- 
mal area dA in region 2 ( "point charge" ) , one obtains 



dF ela = 



dA 



2 Anee 



cr 2 f°° cr 2 

27r/ exp(— nr)dr — dA 

rres J 4se K 



(C3) 



where the assumption of flatness of the surrounding is 
expected to hold everywhere (since we integrate to r = 
oo), but this is of no importance as n effectively cuts 
off the integration. The equation above was obtained 
by summing the screened interaction of a point charge 
(infinitesimal patch) with the rings of infinitesimal width 
and integrating over all the rings in the plane. Note how 
the electrostatic energy of the patch of area dA does not 
depend on its position, which means that the total energy 
of charges in region 2 is 



F 



1:2 



4e£o« 



-A, 



(C4) 



where A 2 is the area of region 2. 

To calculate the electrostatic free energy of region 1 
we approximate it as a cylindrical shell with radius of 
R sin #o . We introduce the coordinate z measured along 
the cylindrical shell, going from z — at the contact 
of regions 1 and 2 to z — 1/k at the edge. The rel- 
evant infinitesimal area is now dA — 2irR sin z9q dz . 



One could introduce a better approximation of region 
1 as a conical shell with the infinitesimal area d^4 = 
2tt(R — zcosi9o)sm'dodz , but this can obviously be 
written as a sum of two parts, the first corresponding 
to the cylindrical shell area above, and a higher order 
contribution which is irrelevant for our purposes. 

The electrostatic energy of region 1 can now be calcu- 
lated as follows: 



F e ll — 



2irRsmd a 2 



Z—l I K 



dz 



-dS z 



(C5) 



where S z now represents the surface obtained by an in- 
tersection of a circle, centered on a certain z coordinate 
and of radius 1/k, and the cylindrical shell. This surface 
is a circle without a circular segment of angle 20o, where 
0o = arccos(l — kz). The integrations over the surface 
can be written as 



-dS z = 2 







/ d 


*r 




p r — oo 


/ 

/0=7r/2 


/ 

Jr=0 



e- Rr dr , (C6) 



where the limits of integration in the second integral 
were extended to infinity, in the spirit of the deriva- 
tion of electrostatic interaction in the region 2. The lim- 
its of integration are thus either terminated at the edge 
l^max — (1 — z k)/(k cos 4>)] or extended to infinity wher- 
ever possible, i.e. where the local neighborhood is similar 
to the one in region 2. The second integral (7a) is easy, 
but the first one (7i) is still nontrivial. We write it down 
in full, where we take into the account the integration 
over z as well: 



h 



1/k 



dz 



tt/2 



r dr. 



(C7) 



One can interchange the order of integration over z and 
(f> , since the integration limits are not interdependent. 
Carrying out the integration over r by introducing l — 
kz = cos <po as a new variable, while defining 



B = I d0i cos < 



l — exp — 



l 



cos < 



we get 



h 



2 /7T 

^2 V2 



-B 



(C8) 



(C9) 



Putting together both integrals in Eq. IC6I we see that 
in comparison with the same integral in region 2, which 
yielded ttk~ 2 , the corresponding contribution in region 
l is reduced by a factor of 2£>/7r because of the effects 
of the edge. The electrostatic contribution to the free 
energy of region l is thus 



27ri?sintf CT 2 A B 
7 e /,i = 5 1 



(CIO) 
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The total electrostatic energy of the incomplete shell Numerical integration gives B w 0.726379. We define 
is then 

F e i = -r^—A 2 + F eL1 , (Cll) / = A » 0.115607 , (C14) 

where A 2 is the area of region 2, 

and thus 

A 2 =27ri? 2 (l + costf ) - — sini? • (C12) 

7ri? 2 cr 2 /l + cosi9 / 



Therefore, we get for the total electrostatic energy ^el — se(jK ( 2 /ci? s * n ^° j ' 



^ei = 7, o 5 Sln1? • (C13) 
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